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\  Abstract 

The  set  of  all  bivariate  positive  quadrant  distributions  is  neither  compact 
nor  convex.  But  the  set  of  all  bivariate  positive  quadrant  distributions  with 
fixed  marginals  is  a  convex  set.  These  convex  sets  are  compact  in  the  case  of 
discrete  bivariate  distributions  if  the  marginals  have  finite  support.  A  simple 
method  to  enumerate  the  extreme  points  of  these  convex  sets  is  given. 

In  the  context  of  contingency  tables  for  testing  the  null  hypothesis  indepen¬ 
dence  against  the  alternative  positive  quadrant  dependence  one  can  use  the  method 
of  extreme  point  analysis  to  compare  the  performance  of  various  tests. 


AMS  1980  subject  classification:  Primary:  62H05  Secondary  62H17 

Key  words  and  phrases:  Extreme  point,  convex  set,  compact  set,  bivariate  distri¬ 
butions,'  positive  quadrant  dependence,  negative  quadrant  dependence,  power  function 
and  contingency  tables. 


1.  Introduction.  Let  be  the  collection  of  all  probability  measures  defined 

2 

on  the  Borel  a-field  of  the  two-dimensional  Euclidean  space,  R  .  For  p  in 

let  p  and  P  denote  the  marginal  probability  measures  on  the  Borel  o-field 

1  2 

of  R  ,  i.e.,  U1(A)  =  p( A  x  R)  and  pg(A)  =  p(R  x  A)  for  every  Borel 

subset  of  R  .  A  p  in  is  said  to  be  positively  (negatively)  quadrant  dependent 

(PQD)  (NQD)  if  ®)  x  (d,  «)}  >_  (<_  )  y^fc,  *))  P2([d,  for  6111  c’  d 

in  R  .  Let  M(PQD)  denote  the  collection  of  all  probability  measures  in  Mg  which 
are  PQD  and  M(NQD)  the  collection  of  those  in  Mg  which  are  NQD  .  The  main 
objective  of  this  paper  is  to  study  the  structure  of  M(PQD)  and  M(NQD)  .  Using 
the  structure  theorem  given  in  Section  3,  we  give  some  applications  in  the  area  of 
testing  hypotheses  in  Contingency  Tables  in  Section  U.  For  further  information  on 
properties  of  measures  which  are  PQD  or  NQD,  see  Lehman  or  Eaton  ^  . 


2.  Preliminary  results .  The  sets  M(PQD)  and  M(NQD)  are  not  convex  sets. 
Examples  are  easy  to  provide.  We  give  a  simple  sufficient  condition  under  which 
convex  combinations  of  measures  in  M(PQD)  are  again  in  M(PQD).  Let  M  denote 
the  collection  of  all  probability  measures  on  the  Borel  a  -field  of  R.  For  v  ,  n 
in  ,  v  is  said  to  be  stochastically  smaller  than  p  if  >{  &  ,")}<_  Tl(  [c  ,  oo 
for  all  c  in  R.  If  v  is  stochastically  smaller  than  r)  ,  we  use  the  notation 
v  <  st  ii  .  The  relation  "  <_  st  "  is  a  partial  order  on  M^.  For  y  ,  X  in  M 
we  say  that  y  is  stochastically  smaller  than  X  if  y^  st  X^  and 

y^  £  st  Xg-  In  this  case,  we  again  use  the  same  notation  y  <  st  X  .  The 

relation  "  <_  st  "  on  is  only  a  pre-order.  See  Definition  2.1  of  Eaton  £  , 

p.ll] . 

Proposition  1.  If  y,x  e  M(PQD)  and  either  y  <_  st  x  or  X  <  st  y  , 

then  ay  +  (l  -  a)x  e  M(PQD)  for  any  Ck  a  <_  1. 

Proof.  We  need  to  show  that  Jay  (l  -  a)TJ{|jc,  »)  x  [d,  »)} 

Ca,Ji  +  (1  -  a)  xjf  [<?,<*>)}  [ay 2  +  (l  -  a)  Xgl'f  )  for  all  c,  d  in  R.  We 

note  that  £ay  +  (1  -  a)^{Q,  ®)  x  [d,  °°)}  >_  00 ^  y2{[d»  00  ^  + 

(1  -  a)  X.^^,  <*>)}  Xg{  [d ,  ®)}  and  ay1{[c,  »)}  y2([d,  *))  + 

(1  -  a)  x1([c,  »)}  x2{[d,  »)}  -  a2  y^Jc,  *))  y2^[^»  *)) 

a(l  -  a)  y ^ { [_c ,  »)}  X2([d,  »)}  -  a(l  -  a)  y2([d,  X^jc,  °°)} 

(1  -  a)2  X1([c,  »)}  X2{|d,  <*>)}  =  a(l  -  a)[^  y-^Jc,  y2{[d»  00 ^  + 

X-^J^c,  «)}  Xg{  [d,  »)}  -  y2([d,  ®}}  X1([c,  «)}  -  y1([c,  °°)}  X2([_d,  »)}J 

=  a(l  -  a)  Cy^  [c ,  ®)}  -  X1([c,  00  ~  X2(Q,  °°)7]  >_  0,  if 

y  <  st  X  or  X  <  st  y  .  This  completes  the  proof. 

An  analogous  result  holds  for  measures  which  are  NQD  as  declared  in  the 
the  following  proposition. 


Proposition  2.  Let  y,X  e  M(NQD).  If  y  x  1  st  X  1  and  X  2  £  st  y  or 
X  ±  <_st  y  ^  and  y  g  <_  st  X  then  ay  +  (l  -  a)  X  e  M(NQD)  for  any 
0  <  a  <  1. 

Corollary  3 .  Let  y  ,X  e  have  the  same  marginals,  i.e.,  y  =  X  and 

U  2  =  X  g.  If  y  ,X  e  M(PQD),  then  ay  +(l-a)Xe  M(PQD)  for  any 

0  <  a  <  1.  If  y,X  e  M(NQD),  then  ay  +  (l  -a  )X  e  M(NQD). 

In  order  to  determine  the  structure  of  M(PQD)  and  M(NQD),  Corollary  3 

provides  a  pointer.  Let  y  ^  and  y  be  two  given  probability  measures  on  the 
Borel  c  -field  of  the  real  line  R.  Let 

MpQD^  y  -p  y  2)  =  (  X  e  M(PQD)  ;  X  i  =  y  i  and  A  2  =  y  2  ^  ’ 

and 

mNqD(  y  ]_» y  2)  =  {  X  e  M(NQD)  ;  X  1  =  p  1  and  X  2  =  w  2  ^  ’ 

Mp(^D( y  ^,y  is  the  collection  of  all  those  probability  measures  X  which  are 
PQD  and  having  the  same  marginals  y  ^  and  y  ^ •  This  set  is  more  tractable  than 
M(PQD)  as  the  following  corollary  points  out. 

Corollary  U>  MpQD^yl’p2^  and  ^NQD ^  y  1 ’ w  21  are  convex  sets. 


Our  next  objective  is  to  determine  the  extreme  points  of  the  above  convex 
sets.  If  Mp^(  y  ^,y  is  a  compact  set  in  some  relevant  nice  topology,  then  one 
can  write  every  X  in  Mp^(  y  ^,y  as  a  mixture  of  its  extreme  points.  See 
Phelps  [  ].  Then  a  knowledge  of  all  extreme  points  of  the  set  Mp^(  y  ^,y  provides 
us  an  insight  into  the  structure  of  Mp^(  y  ^ ,  y  ^ ) .  In  the  sequel,  we  concern 
ourselves  exclusively  with  the  convex  set  y  ,  ,y  .).  The  convex  set  M.T._(  y  ,  ,y 


A 


can  be  studied  in  a  similar  vein.  The  determination  of  all  extreme  points  of  the 
set  ^2^  seems  to  be  difficult  in  the  general  case.  However,  in  Section  3 

we  give  a  complete  description  of  all  extreme  points  of  MpQD(  ^  y2)  when  y1 

and  are  some  special  discrete  measures  with  finite  support.  The  method  given  in\ 
Section  3  can  be  adapted  for  any  discrete  measures  ^  and  ^  with  finite 
support.  Before  sailing  into  the  next  section,  we  identify  one  particular  measure 
as  an  extreme  point  of  Mp^(  Vg)  in  the  general,  case.  For  probability  measures 

^  and  ^  on  the  Borel  a-field  of  the  real  line,  let  V  x  ^g  denote  the 
product  probability  measure  in  Mg. 

Proposition  5.  The  product  probability  measure  h  x  ^g  is  an  extreme  point 

of  MPQD(  V  V  and  0150  of  W  V  V- 

Proof.  One  can  prove  this  result  by  using  the  elementary  fact  that  if  c  1.  a, 
d  ia,  0  <  a  <  1  and  ac  +  (l-a)d  =  a,  then  c  =  a  and  d  =  a,  and  the 

fact  that  any  measure  in  Mg  is  determined  by  its  values  on  sets  of  the  form 
lc,  00 )  x  Id,  °°)  ,  c,  d  real. 

Remark .  The  above  result  is  in  marked  contrast  to  the  following  result  of  Kemp  [3], 
Let  M(  y  ^ ,  Mg)  denote  the  collection  of  all  probability  measures  \  in  Mg 
satisfying  =  y^  and  Xg  =  y  .  Then  M(  y^,  yg)  is  a  convex  set  ,  and 

V1  x  Ug  is  an  extreme  point  of  M(  y  ,  y  )  if  and  only  if  either  y^  is 
0-1  valued  or  Pg  is  0-1  valued. 


3.  Extreme  points.  In  this  section,  we  assume  that  the  probability  measures 


y  and  y  have  finite  support.  For  simplicity,  assume  that  the  support  of  y 
1  2 


is  {1,2,3, ... ,m}  and  that  of 


is  {1,2,3, ... ,n} 


Any  probability 


measure  in  M2  with  support  contained  in  { (i , j )  ;  i  =  1  to  m  and  j  =  1  to  n} 


and  with  marginals  y  and  y  can  be  identified  as  a  matrix  P  =  {p..} 

12  ij 


of  order  m  x  n  satisfying  p„  >_  0  for  all  i  and  j,  ^-2  Pfj  =  P  ({j)) 


for  all  j  =  1  to  n,  and 


2 .  ,  p..  =y  ({i})  for  all  i  =  1  to  m.  Let  p  = 


“j=l 


y  ({ i} ) ,  i  =  1  to  m  and  q.  =  y  ({j}),  j  =  1  to  n.  The  objects  of  interest  are 
1  J  2 


the  convex  sets  Mp^Cy  )  and  )  •  We  denote  M^Cy  )  by 

MpqD<P1,P2.---.Pn,;<Il.<l2 . MNQD(,,l-p2--'-''>„;<1l'q2 . V  '  TheSe 

convex  sets  are,  obviosly,  compact.  Moreover,  they  have  each  a  finite  number 


of  extreme  points.  Each  of  these  sets  is  determined  by  a  finite  number  of  hyper- 
,nri-n 


planes  in  R  .  Now,  we  describe  a  method  of  generating  extreme  points  of 
MPQD^P1,P2* "  '  *  ,Pm’ql’q2’  ’  ‘ '  ,C*n^  dn  some  sPecial  cases  of  m  and  n.  The  method 


so  described  generalizes  to  for  any  m  and  n. 


The  case  of  2x2  tables . 

Let  p,  q  be  two  numbers  in  the  open  interval  (0,  1).  Let  p^  =  p, 

P2  =  1  -  P,  =  q  and  q£  =  1  -  q.  A  matrix  P  =  {p^}  °f  order  2x2  with 

the  properties  p  >_  0  for  all  i  and  j,  p^  +  p12  =  P,  P21  +  P22  =  1  ~  P> 

P11  +  P21  =  q  and  P12  +  p22  =  1  "  q  ls  a  member  of  MpQD^pl’p2’ql’q2)  if 
and  only  if  p 2<1 2  1.  P22  —  P2^^2’  a  ^b  denotes  the  minimum  of  the  two  numbers 

a  and  b.  The  extreme  points  of  this  convex  set  MpQ2>^P2’P2’ql’q2^  are  8iven  by 


*tji  *  v sr^  ir*  jw  aimvv.v  V  VT"LW  < *V=". -  .’777777^^.'^  -^/ 


plql  Plq2 


P2ql  P2q2 


Plql  Plq2 


qx  q2  -  P2Aq2 
0  P2Aq2 


when  p2a  q2  =  P2  5 


P2ql  p2q2 


p2  ~  p2/\q2  p2Aq2 


when  A  q2  =  q2 • 


If  p  =  q  =  1/2,  then  the  extreme  points  of  MpQD(l/2,l/2;l/2,l/2)  are  given  by 


1/A  1/A 


1/A  1/A 


1/2  0 


0  1/2  . 


On  the  other  hand,  the  extreme  points  of  the  set  !:'i/2,l/2;l/2,l/2) ,  the  set  of 
all  matrices  P  =  {p„}  of  order  2x2  with  the  properties  p„  >_  0  for  all 
i  and  j,  pu  +  p12  =  p21  +  p22  =  pn  +  p21  =  p12  +  p22  =  1/2,  are  given  by 


1/2  0 


0  1/2 


The  case  of  2x3  tables 


0  1/2 


1/2  0 


Let  p^,p2;q^ ,q2,q^  be  five  positive  numbers  satisfying  p^  +  p2  = 
q^  +  q2  +  q^  =  1 •  We  now  discuss  the  problem  of  identifying  the  extreme  points 
°f  MpQ^(p^,p2;q^,q2,q^)  •  This  convex  set  is  precisely  the  collection  of  all  ma'ric 


7 


satisfying  pn  +  p12  +  P13  =  Pr  P21  +  p22  +  p23  =  ?2>  ?u  +  p21  = 

p12  +  p22  =  q2’  p13  +  P23  =  q3’  P23  p2q3’  P22  +  p23  P2q2  +  P2q3  and 
each  p  is  non-negative. The  problem  of  identifying  the  extreme  points  of 

this  convex  set  can  be  solved  both  graphically  and  algebraically.  The  following 

results  are  helpful  towards  this  goal. 


Theorem  6.  Let  P  =  (p  )  be  a  matrix  of  order  2x3.  If  P  is  a  member  of 
MPQD(pl’P2;ql’q2’q3)’  th6n 

(D  ?2q3  -  p23  -  p2  Aq3  ’ 

(2)  (p2q2  +  p2q3)  V  P23  1  P22  +  P23  £  P2  A (q2  +  P23)  J 


(3)  P12  q2  p22 ’  p13  q3  “  p23’  P21  p2  p22  “  p23  and 

P11  =  P1  p12  "  p13  =  ql  P21  ‘ 


Proof .  Since  P  is  positively  quadrant  dependent,  p23  >_  P2^3  and  p22  +  p23  > 
P2q2  +  P2q3‘  Since  P  has  the  prescribed  marginals  and  non-negative  entries, 

P23-P2*q3’  P22  +  p23  —  P23’  p22  +  p23  ^  P2  and  p22  +  P23  -  q2  +  P23'  This 

proves  (1)  and  (2).  (3)  is  obvious.  (Note .  For  any  two  real  numbers  a  and 

b,  aVb  means  the  maximum  of  a  and  b.) 


In  the  above  theorem,  the  inequalities  (1)  and  (2)  are  very 
important.  They  help  to  reduce  the  problem  of  identifying  the  extreme  points  of 
MPQD^Pl,P2’ql,q2’q3^  to  3  two-diraensional  graphical  problem.  The  following 
result  amplifies  this  remark. 


Theorem  7.  Let  p22  and  p23  be  two  real  numbers  satisfying 


d)  p2q3 1  p23  i  p2ai3 


(2) 


PnoN/CPa^n  £  Poo  PnO  £  Pt  ^  P  T  O  ^  * 


Let  P12  “  q2  P22’  p13 


=  q3  '  P23  and  P21  P2  “  P22  "  P23’ 


P1  "  p12  p13  ql  “  P21 

°f  MPQD(pl,P2;ql’q2,q3)* 


Pll*  sa^*  and  t*ie  matr:*-x  p  =  (Pj_j)  is  a  member 


Proof.  First,  we  show  that  p1  -  p12  -  p13  =  q1  -  P21*  Note  that  p12  +  p^3  - 

^2  +  q3  ~  ^p22  +  P23^  and  ql  +  P12  +  P13  ^  ^p22  +  P23^  P1  +  P2  +  P21  ~  P2 

=  pl  +  p21.  Consequently,  p.^  -  (p^2  +  P13>  =  -  P21»  Next,  we  show  that 

each  p^  is  non-negative.  From  (1),  P23  >  0  and  q^  -  p23  =  p^3  0 .  From  (2)J 

p22  -0'  P2  ”  p22  "  p23  =  P21  -  °  and  q2  ‘  p22  =  p12  -  °‘  Xt  remains  to  be 
shown  that  p^  is  non-negative.  It  is  clear  that  P11  +  P12  +  p13  +  P21  +  P22  + 

P23  =  1.  It  suffices  to  show  that  p^2  +  p13  +  p21  +  p22  +  p23  <_  1.  Note  that 

p12  +  pi3  +  P21  +  p22  +  p23  =  q2  +  q3  "  (p22  +  p23}  +  P2  "  (p22  +  p23}  +  (p22  +  P23] 

=  q2  +  q3  +  P2  -  (P22  +  P23)  1  q2  +  q3  +  P2  '  P2q2  ”  p2q3  (by  (2))  =  q2  +  q3 

+  P2(l  -  q2  -  q3)  -  1  -  q1  +  P2qx  =  1  "  q^l  ~  P2)  -  1  "  P^  <  1-  ^  shows 

that  p^j  is  always  strictly  positive. 

It  is  obvious  that  the  matrix  P  =  (p^)  has  the  prescribed  marginals 
P;L,P2  and  q;L,q2,q3.  From  (1)  and  (2),  it  follows  that  P  is  positively 
quadrant  dependent . 


The  purport  of  the  above  theorem  is  that  it  is  enough  to  determine  two 
numbers  p22  and  p23  satisfying  (1)  and  (2)  from  which  we  can  build  a 
matrix  P  =  (p±  )  belonging  to  MpQ^P^P^i’^’^  *  We  cal1  (1)  and  (2) 
to  be  the  determining  inequalities.  This  theorem  also  provides  modus  operandi 

to  construct  extreme  points  of  ^pQD^pi,p2’ql,q2’q3^  ‘  Intu*t*ve-*-y »  one  can 
proceed  by  looking  at  those  pairs  p22  and  p23  of  numbers  for  which  equality 
holds  in  (1)  and  (2)  either  on  the  right  or  on  the  left  to  yield  extreme 
points.  Tl^s  method  would  not  exhaust  all  the  extreme  points  of  ^pqd^pi ,p2 ’ql ,q2 ,q3* 


The  reason  is  that  the  bounds  in  (2)  are  variable  and  depend  on  the  value  of 


we  can  Partition  the  interval  tP2q3 *  p2Ac^  into  sub-intervals  so  that 
the  bounds  in  (2)  become  stable  over  each  sub-interval,  then  we  can  isolate 
the  extreme  points.  We  explain  this  method  graphically  and  algebraically  as 
follows . 

Graphical  solution 

Determining  inequalities 

(1)  p2q3  -P23  -P2Aq3  * 

(2)  P23V(P2q2  +  p2q3)  -p22  +  P23  lp2A(q2  +  p23^  ‘ 

Case  1  p2q2  +  p2q^  L  P2A  ^3  and  q2  Lp2-  Then  the  determining  inequalities 

become 

(1)  p2q3  i  p23  !  p2A  q3  and 

(2)  p2q2  +  p2q3  1?22  4-  p23  ip2. 

These  inequalities  yield  the  following  graph. 
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Comments .  The  extreme  points,  in  this  case,  are  obtained  by  setting  each 
expression  in  the  centre  of  (1)  and  (2)  equal  to  the  quantities  either  on 
the  left  or  on  the  right  of  (1)  and  (2).  This  will  give  us  4  extreme  points. 

Case  2.  p2q2  +  p2q3  >_  P2Aq3  and  q2  <  p2- 

Then  either  (A)  p2  -  q2  <_  p2q3  or  (B)  p2q3  <  p2  ~  ^2  *  p2^q3  or 
(C)  p2  ~  q2  —  p2^q3  holds.  The  determining  inequalities  and  the  corresponding 
graphs  then  respectively  work  out  as  follows. 


Case  2 (A) .  Inequalities 


U)  p2q3  -  p23  -  P2  Aq3  and 

(2)  p2q2  +  p2q3  £  p22  +  p23  <  P2  . 

Graph.  Same  as  in  Case  1. 


Case  2 (B) .  Inequalities 

(1)  P2q3  -  p23  ~  p2"  q2  and 

(2)  p2,2  +  p2q3  <  p22  +  P23  <  q2  +  P23 

or, 

(1) '  p2  -  q2  <  P23  <.  P2  A  ^3  and 


Comments .  can  be  above  ?2^2  +  ^2^3*  num^er  °f  extreme  points 

is  four. 

Case  3  P2q2  +  P2<J3  —  p2^  ^3  anc*  p2  —  q2’ 

The  determining  inequalities  become 


P2q3  P2q2+P2q3  p2A  q3  p2  P23 

Comment .  The  number  of  extreme  points  is  five. 


Case  4 


P2q2  +  P2q3  <  P2^  q3  and  P2  >  q2’ 

Then  either  (A)  ?2  -  q2  >_P2/\q3  °r(^2q3  <  P2  -  q2  <  P2A  q3 

P2  -  3 2  —  p2q3  holds*  We  now  discuss  each  of  these  cases. 


or 


Case  4(A) .  Inequalities 

(1)  P2^3  i.  P23  —  P2^2  +  P2q3  and 

(2)  P2<l2  +  P2,3SP;,2  +  P231q2  +  p23  , 

or, 


(C) 


jrr 


P2q2+P2q3 


p2  Aq3 


P2q3 


P 


r 

23 


<^2  can  be  below  P2^2  +  p2q3‘  There  are  ^ive  extreme  points. 


Assume,  without  loss  of  generality,  that  P2  ”  q2  —  P2^2  +  P2q3* 

* 

Inequalities 

(!)  p2q3  i  p23  i  p2  -  q2  and 

(2)  P2q2  +  p2q3  l  p22  +  p23  iq2  +  p23  , 

or, 

(1)  P2  -  q2  <  P23  1  P2q2  +  p2q3  and 

(2)  p2i2  +  p2q3  1p22  +  p23  ip2  , 

5 


p2q2  +  **2«*3  '  p23  -p2A<>3 


and 


>.  Jtm  . 


-M 


i.  m ' M  » 


Algebraic  solution 


The  expression  (p2q2  +  P2q3^  ^23  inc*uces  a  partition  of  the 
interval  [p2q3,  P2A  Q3I  into  two  sub-intervals.  These  intervals  are  given  by 

A1  =  ^p23  e  tp2q3*  p2Aq3^  5  <p2q2  +  p2q3^Vp23  =  P2q2  +  p2q3^ 


^2  =  ^p23  ®  tp2q3»  p2 ^  ^3 1  »  (p2q2  p2q3^^P23  p23^  * 

A^  and  A^  may  have  a  common  end  point.  In  a  similar  way,  the  expression 
P2A(q2  +  P23>  induces  a  partition  B3  and  B2  of  the  interval  [p2q3>  p2A  Q-jl 
into  two  sub-intervals.  The  partitions  { A^ ,  A2>  and  {B^,  B2>  of  [p2q3>  p2Aq3] 
induce  a  finer  partition  {C^,  C2,  C3>  of  [p2q3>  p2Aq3l-  Let 


C1  "  ^p2q3 *  al^  ’  C 


[a1,  a2]  ;  C,  =  [a9,  p,AqJ  . 


2>  ^2  43j 


We  now  consider  three  sets  of  inequalities. 


(A)  (1)  p2q3  <_  p23  <_  a3  and 

(2)  P23  V  Cp2q2  +  p2q3^  -  p22  +  p23  -  p2 A  (q2  +  P235 ' 

(B)  (1)  a3  l.P23  —*2  and 

(2)  p23  V (P2q2  +  P2q3^  -  p22  +  p23  -  p2  A  ^q2  +  P23) ’ 

(C)  (1)  a2  i  P23  iP2Aq3  and 

(2)  P23V(p2q2  +  P2q3)  1P22  +  P23  lP2A(q2  +  P23)  . 

In  each  set,  in  the  inequality  (2),  the  max  and  min  symbols  V  ,  A  would 
disappear  after  due  simplification.  From  each  set,  the  central  expression 
in  each  inequality  is  set  equal  to  the  quantity  either  on  the  right  or  on  the 
left.  These  equations  are  solved  for  unknown  values  of  p22  and  P23*  ^he  set 
of  matrices  so  obtained  after  using  Theorem  7  contains  the  set  of  all  extreme 
points  of  MPQD<P1»P2i<l1»cI2»<:l3>  * 


<2)  W  <p2Vl  +  P2V  i  p2„-l  +  P2„  1  P2  A  <Vl  +  p2n>  ' 

<3)  <p2n-X  +  ^  V  <P2V2  +  p2Vl  +  P^n*  -  p2n-2  +  p2n-l  +  p2„  - 

P2  A  C  V2  +  P2„-l  +  Pxn*  ’ 


(n-1)  (p23  +  p24  +  ...  +  p,T,)V(p9q9  +  P9q9  +  ...  +  P9qri)  1 


2n  T  ^2M2  '  K2M3 


P22  +  P23  +  •••  +  P2n  iP2A(12  +  p23  +  P24  +  "•  +  p2l,> 

These  inequalities  provide  a  basis  for  generating  members  of 
JiPQD^l,^2;^l*t*2’  ‘  *  *  ,qn^  ‘  ^°^^ow^n8  theorem  provides  a  method  for  finding 

members  of  mPqD(P1»P2^i><52’ “  *  ,<5n^  * 

Theorem  8  Let  p22,  p23>  ...  ,  p2n  be  (n-1)  numbers  satisfying  the 
inequalities  (1),  (2),  ...  ,  (n-1)  above.  Let  p^  =  q_.  -  P2j »  J  =  2,3,... 
and  p21  -  p2  -  (p22  +  P23  +  ...  +  p2n)-  Then 


’i  -  p2i  '  pi  -  (pi2  +  p13  +  •••  +  pln>  -  pn.  »®y. 


and  the 


matrix  P  =  (p  )  of  order  2  x  n  belongs  to  MpqD(P1  ,P2 ^ ,q2 » •  ■  • , 


Proof .  It  is  similar  to  the  one  presented  for  Theorem  7. 


Now,  we  set  about  obtaining  extreme  points  of  Mp^D(p3,p2;q3>q2, . . . 


algebraically , 
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We  develop  some  sets  of  inequalities  equivalent  to  the  inequalities 
(1) ,  (2) , . . . , (n  -  1)  described  above  with  the  following  properties. 

(i)  Each  set  contains  (n  -  1)  equalities 

(ii)  The  max  and  min  symbols  V  , A  disappear  from  each 
side  of  each  inequality  in  each  set. 

Step  1.  The  interval  [ P 2clri »  AqR]  variation  for  p^n  is  split  into 

3  sub- intervals  using  the  inequality  (2)  .  The  expression 
P2n>/(P2qn_1)  +  P2qn>  Sives  a  Partition  of  the  interval  tp2cln»  P2Aqn3 
into  two  sub-intervals  as  follows.  Let 

A1  *  {p2n  e  lp2qn>  p2Aqnh 

p2n  V*p2Vl  +  p2qn>  ‘  p2n> 

and 

A2  -  {p2n  £  Ip2qn.  P2Aqnl; 

p2n  V(P2’n-l  +  p2q„)  *  p2qn-l  +  p2qn> ' 

Note  that  and  are  closed  intervals  with  possibly  a  common 
end  point  and 

AlUA2  =  tp2V  *2  Hi1* 

(One  of  the  A^'s  could  be  a  null  set.) 

Similarly,  the  expression  p2A(<?n  ^  +  P2n)  8^ves  3  partition  {B^,B2) 

of  the  interval  [P2qn>  ?2^qn^  ^nt0  two  cl°sed  sub-intervals.  The 

partitions  {A^,A2)  and  {B^,B2)  induce  a  finer  partition  {C^C^C^}  of 

[p„q  ,  p„/\  q  1  such  that  each  C.  is  either  a  closed  interval  or  a  null 
r2  n  2'  n  l 

set  and 

ClUC2UC3  '  |P2V  P2/'qnl- 

For  example,  if  A^  =  [a^,a23,  A^a^a^],  =  [b^,b2],  B2  =  [b^b^]  with 

a.  =  b.  =  p„q  ,  a_  =  b_  =  p~A  q  and  a„  <  b„,  then  we  can  take  C.  =  [a.,,a2]. 


*“2  ~  la2’^2^  anc*  ^3  =  [^2*^3  *  ma^ntain  a  conformable  notation, 

write  ax  =  a^,  &2  =  aj.^,  &2  =  a[2],  b2  =  a'[2],  b2  =  a[3^  and  b3  -  a 
we  have 

[p2qn’  p2A^n3  =  [a[l]’  a[l]]  [a[2],ai2]](J[a[3],a[3]] 

where  a^  =  p2qn>  aj^  =  a[2],  aj2]  =  a[3]  and  aj3]  =  P2/1  V 
The  inequalities 

(1)  p2qn  -  p2„  i  p2qn  and 

(2)  f2lV(p,Vl  +  P2qn)  <  p2n.1  +  p2n 


i  p2A<Vl  +  p2n> 

are  equivalent  to 


[1] 

a[l]  ip2nia[l]  and 

[13(2) 

p2nV(p2qn-l  +  P2V  i  p2n-l  +  p2n 

i  P2  A  <Vl  + 

or 

[2] 

at2]  i  p2n  i  a[2]  and 

[2] (2) 

p2nV<p2Vl  +  p2qn)  i  p2n-l  +  p2n 

i  p2^  (Vl  +  p2n^  ’ 


or 

[3] 

a[3]  i  P2„  i  a[3J  and 

[3] (2) 

p2n'/(p2Vl  +  P2qn)  ±  P 

1  P;A(Vl  +  p2n)  ‘ 


In  each  of  the  inequalities  [1](2),  [2]  (2)  and  [3]  (2),  the  symbols 

V  and  A  would  disappear  when  p2^  is  confined  to  the  limits  imposed  by 
[1],  [2]  and  [3]  respectively.  To  mark  this  freedom  from  the  symbols 

V  and /j  ,  let  us  rewrite  the  above  inequalities  as  follows. 
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[1] 

a[l]  -  P2n  -  a[l] 

and 

[1] (2) 

a[l] (2)  -  p2n-l  + 

P2n  -  a[l] (2) » 

or. 

[2] 

a[2]  -  P2n  -  a[2] 

and 

[2] (2) 

a[2] (2)  —  P2n-1  + 

P2n  -  a[2] (2) ’ 

or. 

[3] 

a[3]  —  P2n  —  a[3] 

and 

[3] (3) 

a[3]  (2)  -  P2n-1  + 

P2n  -  af 3] (2) * 

Step  2.  Using  the  inequality  (3),  we  sub-divide  each  of  the  intervals 
given  by  the  inequalities  [1](2),  [2]  (2)  and  [3] (2)  onto  three  sub¬ 
intervals  or  inequalities.  Let  us  concentrate  on  [1](2).  The  procedure 
outlined  below  is  similar  to  the  one  given  in  Step  1.  The  expression 

(p2n-l  +  p2n)v'(>’2V2  +  p2qn-l  +  p2qn>  bl£urcates  the 
inequality  [ 1 ] (2)  into  two  inequalities.  The  sets 


D1  '  (p2n-l 

+  p2n  « 

:  [atl] 

(2)'  a[l] (2) ] ; 

(p2n-l 

+  P2n> 

v'(P2qn-2  +  p2qn-l 

+  P2q„) 

=  P2n- 

1  +  p2n}  and 

D2  "  tp2n-l 

+  p2n  ' 

£  [a[l] 

(2)’  a'[l](2)]; 

(p2n-l 

+  p2n^ 

V(p2qn-2  +  p2qn-l 

+  p2q„> 

■  p2qn-2  +  p2qn-l  +  ^n1 

do  give  the  necessary  bifurcation  of  the  inequality  [ 1 ] ( 2 )  with  P2n 
satisfying  the  inequality  [1],  Similarly,  the  expression  ^n_2  +  p2n-l  +  p2n^ 
bifurcates  the  inequality  [ 1 ]  (2 )  into  two  inequalities.  These  two  bi¬ 
furcations  would  give  three  inequalities  equivalent  to  [ 1  ]  ( 2 ) .  Let  these 
inequalities  be  denoted  by  [11],  [12],  [13].  Let  us  write  these  inequalities 


including  [1]  as  follows. 


[1J  a(1]  <  p2n  <  ajjj.  _ 

til]  «[u]  <  P2n-i  +  J>2„  i  *iu]' 

or 

[U]  a[i2]  —  P2n-!  +  P2n  —  a|12]» 

or 

1131  a  [13  J  —  p2n-l  +  p2n  —  a[l3]’ 

where  a|^j  ~  a[i2]  an£*  a[l2]  =  a[13]‘  bifurcating  the  inequality  [1] 
further,  necessary,  ensure  that  a'^j  <  a[ll]  aS  p2n  —  p2n  1  +  p2n  *s  t0 
be  guaranteed.  This  procedure  is  applied  to  each  of  the  inequalities 
[2]  (2)  and  [3]  (2). 

We  can  now  claim  that  the  inequalities  (1),  (2)  and  (3)  are 
equivalent  to  the  following  9  sets  of  inequalities.  (We  describe  these 
inequalities  using  a  suggestive  notation.) 


[1] 

“ill  i  p2n  -  a[X]’ 

[11] 

a[ll]  i  p2n-l  +  p2n  i  aill]>  and 

[11] (3) 

Cp2n-1  +  p2n> '' <p2qn-2  +  p2qn-l  + 

p2qn> 

i  p2n-2  +  p2n-l  +  p2n 

1  p2A  <Sn_2  +  P2n-1  +  p2n>  ’ 

or 

[1] 

a[l]  ip2nia(l]> 

[12] 

a[12]  -  p2n-l  +  p2n  -  a[12]  and 

[12] (3) 

(p2i>-1  +  p2n)'/(p2V2  +  p2qn-l  + 

p2q„> 

1  p2n-2  +  p2n-l  +  P2n 
±M(qn-2+P2n-l+P2n)’ 


[13]  (3)  (P2n_i  +  P2n),/(p2q 
1  p2n-2  +  P2: 

i  p2  a (qn_2  ' 


or. 

[2] 

a[2]  -  p2n  -  a[2]  ’ 

[21] 

a [ 21 ]  -  p2n-l  +  p2n 

[21] (3) 

(p2n-l  +  p2n)v,(p2qi 

1  p2n-2  +  p2i 

ip2A(qn-2  ‘ 

* 

u 

o 

[2] 

a[2]  -  p2n  -  a'  [2]  ’ 

[22] 

a [ 22 ]  -  P2n-1  +  P2n 

[22] (3) 

(p2n-l  +  p2n>  V  (p2q, 

i  P2n-2  +  P2i 

i  p2  A(qn-2  ' 

or. 

[2] 

a[2]  -  P2n  -  a' [2]  ’ 

[23] 

a [ 23 ]  -  P2n-1  +  P2n 

[23] (3)  (p. 


+  PoJ  ^(p9q 
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or. 


[3] 

aI3]  —  p2n  —  a[3] 

[31] 

al31]  —  p2n-l  +  p2nia’[31]  and 

[31] (3) 

(p2n-l  +  p2n>  '"<p2<>n-2  +  p2  Vl  +  P2%,> 

i  p2n-2  +  p2n-l  +  p2n 

ip2A(V2  +  p2n-l  +  p2n^  ’ 

or. 

[3] 

a[3)  -  p2n  -  a[3] 

[32] 

a[32]  i  p2n-l  +  p2n  i  a[32]  and 

[32] (3) 

(p2n-l  +  p2n)  '/<P2V2  +  p2Vl  + 

i  p2n-2  +  p2n-l  +  p2n 

i'2<l(V2  +  p2n-l  +  p2n^  ’ 

or. 

[3] 

a[3]  —  p2n  —  a[3] 

[33] 

a[33]  i  p2„-l  +  p2n  i  ai33)  and 

[33] (3) 

<p2n-l  +  p2n>  V  (p2V2  +  p2\.-l  +  P2V 

i  p2n-2  +  p2n-l  +  p2n 

ip2''(V2  +  ',2Vl+P2n>' 


The  following  facts  should  be  borne  in  mind  before  proceeding  to 
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(iii)  In  each  of  the  nine  inequalities  [11]  (3),  [12] (3), 
[133(3),  [21] (3) ,  [22] (3) ,  [23]<3),  [31] (3) ,  [32] (3) 
and  [33] (3),  the  symbols  V,A  would  disappear  as  long 
as  P2n  and  P2n_i  +  P2n  are  confi-ned  to  the  limits  of 
the  corresponding  inequalities  within  each  set. 

(iv)  The  inequalities  [1],  [2]  and  [3]  obtained  in  Step  1 
can  further  be  bifurcated  to  accomodate  (iii) .  This  is 
helpful  in  conjunction  with  the  inequalities  [11],  [12], 
[13],  [21],  [22],  [23],  [31],  [32]  and  [33]  in  each 
appropriate  category  to  get  rid  of  the  symbols  V  and  A 
in  the  inequality  (3). 

Step  3.  Using  the  inequality  (4),  we  split  the  last  inequality  in  each 
set  above  into  three  inequalities.  The  procedure  is  similar  to  the  one 
explained  in  Step  2. 


Step  (n  -  2)  .  Using  the  inequality  (n-1),  we  split  the  last  inequality 
in  each  set  into  three  inequalities.  At  the  beginning  of  this  step, 
each  set  has  (n  -  2)  inequalities. 

At  the  conclusion  of  Step  (n-2),  the  inequalities  can  be  arranged 
pyramidically  as  in  the  following  diagram,  assuming  that  at  every 
stage  we  split  each  inequality  into  three  inequalities. 
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Each  inequality  at  the  final  Step  (n-2)  is  coupled  with  the  in¬ 
equality  (n  -  1) .  Thus  we  have  3°  2  sets  of  inequalities  with  each  set 
containing  (n  -  1)  inequalities  including  the  inequality  (n-1).  In  the 
last  inequality  (n-1),  the  symbols  V  and  /\  would  disappear.  The 
following  are  examples  of  two  sets  of  inequalities  when  n  =  6,  where  at 
every  stage  each  inequality  gave  rise  to  exactly  three  inequalities. 

Example  1. 

W  a[i]  1  P26  -  a[l]» 

[i2]  a[12]  -  P25  +  P26  -  a[12] » 

tl21]  a[121]  -  P24  +  P25  +  P26  -  a[l21]* 

[1213]  a[1213]  1  P23  +  P24  +  P25  +  P26  1  a [ 1213 ] 

and 

(5)  (p23  +  P2A  +  P25  +  P26)  V/Cp2q2  +  P2q3  +  p2q4  +  p2q5  +  p2q6) 

1  p22  +  P23  +  P24  +  P25  +  P26 
<  p2/\(q2  +  p23  +  p24  +  p25  +  p26). 

Example  2 . 


[3] 

a[3]  -  P26  -  a[3]  ’ 

[32] 

a [32]  -  P25  +  p26  - 

a [ 32 ] * 

[321] 

a[3211]  -  P24  +  p25 

+  P26  -  3 [3211] * 

[3212] 

and 

a [ 3212 ]  -  p23  +  P24 

+  p25  +  P26  —  3 [ 3212 ] 

(5)  (p23  +  p24  +  p25  +  P26>  ^P2q2  +  p2q3  +  P2q4  +  P2q5  +  P2q6} 

-  p22  +  p23  +  P24  +  P25  +  P26 
1  P2  A(q2  +  p23  +  P24  +  p25  +  P26}- 


Enumeration  of  extreme  points  of  MrqpiPpP;?  q1»q2* • • • »qn) 


Let  denote  the  collection  of  all  those  matrices  P  =  (p^)  in 
MpQD(Pi»P2»  *  *  *  ,qn^  ^or  w^ic^  equality  holds  either  on  the  left 

or  on  the  right  in  all  the  inequality  expressions  in  some  set  of  the 
3n  2  sets  described  above.  Let  Ji#  i  =  2,3,..., n,  denote  the  collection 
of  all  matrices  P  =  (p^)  in  MpQD(P^>P2*  qi,q2,q3^  for  which  strict 
inequality  holds  in  (i  -  1)  inequality  expressions  and  equality  holds 
either  on  the  left  or  on  the  right  in  the  remaining  (n  -  1)  -  (i  -  1)  =  n  -  i 
inequality  expressions  in  some  set  of  the  3°  2  sets  described  above. 


Obviously, 


u 

MPQD^p1,P2’  qi»q2’*’*qn^  =  ^  Ji * 


Theorem  9.  The  set  of  all  extreme  points  of  MpQD(P^*P2»  q1>q2»*,,qn) 


is  a  subset  of  J^. 

Proof.  It  suffices  to  show  that  every  matrix  in  mPqd(Pi»P2»  q1»q2»***qn) 
can  be  expressed  as  a  convex  combination  of  members  of  J^.  Every  member 
of  J^,  trivially,  has  the  above  property.  Let  P  =  (p^j )  be  any  member 
of  J2*  We  express  P  as  a  convex  combination  of  two  members  of  J^. 

Since  P  e  J-,  strict  inequality  holds  in  exactly  one  inequality  expression 
of  some  set  of  inequalities.  For  simplicity,  let  the  strict  inequality 


occur  in  [23] ,  i .e. , 


a  [  2  3  ]  < 


p2n-l  +  P2n  *  a[: 


and  equality  holds  in  the  remaining  (n  -  2)  inequality  expressions  either 
on  the  left  or  on  the  right.  Let  Q  =  (q^)  be  the  matrix  in 
MpQD(pi>p2;  q1»q2»‘**»qn)  f°r  which  equality  holds  on  the  left  in  [23] 
and  equality  holds  in  the  remaining  (n  -  2)  inequality  expressions  in 
exactly  the  same  place  as  that  of  P  in  the  same  set  of  inequalities. 
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Let  S  =  (s t j )  be  the  matrix  in  MpQD(P1»P25  qi»q2*  *  ‘  *  ,qn‘>  for  which 

equality  holds  on  the  right  in  [23]  and  equality  holds  in  the  remaining 

(n  -  2)  inequality  expressions  in  exactly  the  same  place  as  that  of  P  in 

the  same  set  of  inequalities.  The  existence  of  Q  and  S  is  guaranteed 

by  Theorem  8.  We  also  observe  that  q^n  ^  <  P2n_^  <  s2n_i*  There  exists 

0  <  a  <  1  such  that  p„  .  =  aq„  ,  +  l-ct)s«  ..  Indeed,  after  laborious 

2n-l  2n-l  2n-l 

calculations,  one  can  check  that  P  =  aQ  +  (l-a)S. 

Let  P  be  a  member  of  J^.  We  can  express  P  as  a  convex  combination 
of  two  members  of  J^.  This  then  would  imply  that  P  is  a  convex  com¬ 
bination  of  members  of  J^.  Since  P  e  J^,  strict  inequality  holds  in  2 
inequality  expressions  and  equality  holds  either  on  the  left  or  on  the 
right  in  ea^h  of  the  remaining  (n  -  3)  inequality  expressions  of  some 
set  of  inequalities.  For  simplicity,  let  the  strict  inequality  hold 


in  [123]  and  [23112],  i.e., 


[123] 


<  p2n-2  +  p2n-l  +  P2n  "  a[ 


[123] 


a[23112]  <  P2n-4  +  P2n-3  +  P2n-2  +  P2n-1  +  P2n  <  a[23112] * 

Let  Q  «  (q^)  be  the  “atrix  in  mPqDCp1,P2;  qx*q2 . qn)  such  that  q2n  =  P2n’ 

q2n-i  =  p2n-l’  q2n-2  =  p2n-2’  ecJuality  h°lds  on  left  of  [23112]  and  equality 
holds  in  the  remaining  inequality  expressions  in  exactly  the  same  place 
as  that  of  P  in  the  same  set  of  inequalities.  Let  S  =  (s^)  be  the  matrix 
in  MpQD(PrP2;  qrq2 . qn)  such  that  s2n  =  p^,  s^  =  P2n_1# 

s2n_2  *  p2n  2,  equality  holds  on  the  right  of  [23112]  and  equality  holds 
in  all  the  remaining  inequality  expressions  in  exactly  the  same  place  as 
that  of  P  in  the  same  set  of  inequalities.  One  can  check  that  P  =  aQ  +  (l-a)S 
for  some  0  <  a  <  1 . 


%  *»  ‘•  fc*'1  » *»  • 
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Proceeding  this  way,  one  can  show  that  every  member  of  is  a 
convex  combination  of  two  members  of  ^  for  i  =  2,3,..., n. 

This  completes  the  proof. 

The  above  theorem  helps  to  identify  the  extreme  points  of 
MPQD^P1*P2*  qi,q2* '  "  ,qn^  as  follows* 

Step  1.  Create  3n  ^  sets  of  inequalities  as  described  above. 

Step  2.  For  each  set  of  inequalities,  set  the  central  expression  in 
each  inequality  equal  to  the  quantity  either  on  the  left  or  on  the  right. 

Solve  the  resultant  equations  for  P2n»  P2n-1 . P23  an<*  P22* 

Step  3.  Using  Theorem  8,  build  a  matrix  P  =  (p^)  in  MPQD(P1>P2i  *q2 *  *  *  *  * qn^  * 

By  Theorem  9,  every  extreme  point  of  Mpqjj(Pi»P2»  q^.q^* • • • »qn)  arises 
this  way.  There  might  be  a  duplication  of  matrices  on  a  moderate  scale. 

Also,  many  non-extreme  points  might  arise  this  way.  We  can  weed  them 
from  J^.  If  it  is  a  hard  task,  we  can  keep  them.  Retaining  them  will  not 
hamper  us  in  discussing  objectively  merits  of  various  tests  for  contingency 
tables  in  the  next  section.  Once  the  sets  of  inequalities  are  created, 
solving  equations  can  be  programmed  easily. 

Intuitively,  Mpq^Pj/P;)!  q3»q2> • • • >9n)  Can  v3ewec*  as  3  multi-fac¬ 
eted  diamond  and  the  solutions  obtained  above  give  the  points  where 
edges  meet. 

A  numerical  example 

Let  n  =  4,  p^  =  0.6,  p2=0.4,  q^  =  0.3,  q2  =  0.2,  q3=0.2  and  q^  =  0.3. 

Determining  inequalities 

(D  P2q4iP24iP2Aqv 

(2)  P24V(P2q3  +  P2q4)  £  P23  +  P24  1  P2  ^  (q3  +  p24)  * 

(3)  (P23  +  p24)V(P2q2  +  P2q3  +  P2q4)  £  P22  +  P23  +  P24  £  P2^(q2  +  P23  +  P24) 
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Numerically,  these  inequalities  become 

(1)  0.12  _<  p24  <  0.3, 

(2)  P24  V (0.20)  <  p23  +  p24  <  0.4  A  (0.2  +  p24), 

(3)  (P23  +  P24)  V  (0.28)  <  p22  +  p23  +  p24  <  0.4  A  (0.2  +  p23  +  p24) 

The  inequality  (2)  induces  only  two  inequalities  [1]  and  [2] 
equivalent  to  (1) . 

[1]  0.12  _<  p24  _<  0.2  and 

[2]  0.20  <  p24  <  0.3. 

The  inequality  (2)  in  conjunction  with  the  inequalities  [1]  and  [2] 


becomes 


[1]  0.12  <_  p24  _<  0.2  and 

[1  ]  (2)  0.20  _<  p23  4-  p24  _<  0.2  +  p24. 


[2]  0.2  _<  p24  <  0.3  arid 

[2](2)  P24  1  P23  +  p24  -  °*4* 

The  inequality  (3)  induces  two  inequalities  equivalent  to  [ 1 ] (2)  and  also 
two  inequalities  equivalent  to  [2]  (2).  We  have  altogether  4  sets  of 
inequalities  listed  below. 


I  Set 


[I]  0.12  <  p24  <_  0.2, 

[II]  0.20  <_  p23  +  p24  <_  0.28  and 


(3)  0.28  <_  p22  +  p23  +  p24  £  0.4. 


II  Set 


[1]  0.12  <  p24  <  0.2, 

[12]  0.28  £  p23  +  p24  £  0.2  +  p24  and 


p23  +  p24 « p22  +  p23  +  p24 1  °-4- 


Ill  Set 


[2]  0.2  £  p24  £  0.3, 

[21]  P24  1  P23  +  p24  -  °'28  and 

(3)  (P23  +  P24)V  (°-28)  1  P22  +  p23  +  p24  -  °‘4‘ 

IV  Set 

[2]  0.2<p24<0.3, 

[21]  0.28  £  p22  +  P24  1.  0*^  and 

(3)  (p23  +  p24)  V  (0.28)  £  p22  +  p23  +  p24  £  0.4. 

In  the  III  set,  the  number  on  the  left  hand  side  of  [2]  is  larger  than 
the  number  on  the  left  hand  side  of  [21].  We,  further,  bifurcate  [2]  to 
set  right  this  anomaly.  After  this  bifurcation,  we  relabel  the  inequalities 
and  obtain  the  following  sets. 

Set  I 

[I]  0.12  <  p24  <  0.2, 

[II]  0.20  <_  p23  +  p24  £  0.28  and 

(3)  0.28  <  P22+P23+P24  1  0.4. 

Set  II 

[1]  0.12  <  p24  <  0.2, 

[12]  0.28  <_  p23  +  p24  £  0.2  +  p24  and 

(3)  p23  +  p24  ±  p22  +  p23  +  P24  -  °’4' 

Set  III 

[2]  0.2  £  p24  £  0.28, 

[21]  P24  i  P23  +  P24  1  °-28  and 

[3]  0.28  £  p22  +  p23  +  p24  £  0.4. 


.a  ^  v  -  a. 
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Set  IV 


[2]  0.20  £  p24  £  0.28, 


[22]  0.28  <  p23  +  p24  <  0.4 


and 


(3) 


23 


+  P24  1  P99  +  P9^  +  P9A  1  °*4* 


22 


23 


24 


Set  V 


[3]  0.28  £  p24  £  0.3 


[31]  P24  1  P23  +  p24  -  0,4  and 

(3)  p23  +  p24  <  p22  +  p23  +  p24  <  0.4. 

Set  the  central  expression  in  each  inequality  in  each  set  to  the 
quantity  either  on  the  right  or  on  the  left  and  solve  the  resultant 
equations  for  p22>  p23’  p24’  distinct  matrices  in  of  Theorem  9 


are  listed  below. 


(1) 

0.18 

0.12 

0.12 

0.18  1 

(2) 

1  0.3 

0.0 

0.12 

0.18 

0.12 

1 

0.08 

0.08 

0.12 

^  0.0 

0.2 

0.08 

0.12  ( 

(3) 

|  0.18 

0.2 

0.04 

0.18  1 

(4) 

0.3 

0.08 

0.04 

0.18 

^0.12 

0.0 

0.16 

0.12  l 

0.0 

0.12 

0.16 

0.12 

(5) 

0.18 

0.12 

0.2 

0.1 

(6) 

|  0.3 

0.0 

0.2 

0.1 

0.12 

1 

0.08 

0.0 

0.2 

^  0.0 

0.2 

0.0 

0.2 

(7) 

0.18 

0.2 

0.12 

0.1 

(8> 

0.3 

0.08 

0.12 

0.1 

i  0.12 

0.0 

0.08 

0.2 

i  0.0 

0.12 

0.08 

0.2 

(9) 

0.22 

0.2 

0.0 

0.18  1 

(10) 

0.3 

0.12 

0.0 

0.18 

0.08 

0.0 

0.2 

0.12  J 

0.0 

0.08 

0.2 

0.12 

(ID 

0.3 

0.2 

0.0 

0.1  ' 

(12) 

|  0.18 

0.12 

0.2 

0.1 

0.0 

0.0 

0.2 

0.2 

1  0.12 

0.08 

0.0 

0.2 

k^iviviviviY^fiv^iafciyi^ 
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0.2,  0.2,  0.3).  (4)  is  not  an  extreme  point.  For,  we  can  write 

1  2 

(4)  =  -^(2)  +  -^(10).  (15)  is  not  an  extreme  point.  For,  we  can  write 

1  4 

(15)  =  -jr(ll)  +  y(18)  .  (7)  is  not  an  extreme  point.  For,  we  can  write 

(7)  =  -j(3)  +  -^(13).  (8)  is  not  an  extreme  point.  For,  (8)  =  y(4)  +  y(14). 

The  number  of  extreme  points,  therefore,  of  >^^^(0.6,  0.4;  0.3,  0.2,  0.2,  0.3) 
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Now,  it  is  time  to  discuss  the  general  m  x  n  case.  But  the 
arguments  involved  and  the  notation  used  are  so  complex,  it  is  not 
really  worth  the  effort.  Instead,  we  offer  a  thorough  discussion  of 
3x3  case  from  which  one  can  make  a  successful  analysis  of  any  m  x  n 
case.  The  analysis  of  3 x  3  case  is  not  exactly  in  tune  with  the  2x  n 
case.  A  study  of  2x  n  and  3x3  cases  in  tandem  should  help  to  map  out 
a  successful  strategy  to  analyze  mxn  case  in  general. 

The  case  of  3x3  tables 

Let  p^,  p2,  p^>  q2>  ^3’  six  positive  numbers  such  that 

P1  +  p2  +  p3  =  ql  +  q2  +  q3  =  d ^  a  matr^-x  p  =  (pij)  order  3x3 
belongs  to  Mpqj/Pj^ >P2  ,p3;  ql,q2’q3^’  then  °ne  Can  check  that 


(1) 

p3q3  £  p33  £  p3A  q3. 

(2) 

p33  V  (p2q3  +  p3q3}  i  P23 

+  p33 

-  q3  - 

^  (p2  + 

p33** 

(3) 

p33  V  (p3q2  +  p3q3)  i  P32 

+  P33 

—  P3 

^  (q2  + 

P33), 

(4) 

^p23  +  p32  +  p33^  ^  ^p2q2 

+  P2q3 

+  P3q2  + 

P3q3} 

—  P22 

+  p23  + 

■  P32  +  P 

33 

1  (P2 

+  p32  4 

'  p33}  A 

(q2  +  P23  +  P33^ 

hold. 

P  is  PQD  yields  (1)  p^^  : 

-  P3q3 ’ 

(2)  P23 

4  p33  >_  p2q3  4  p3q3. 

(3)  p 

32  +  p33  ^  p3q2  +  p3q3  and 

(4)  p22 

+  P23  + 

P32  +  p33  -  P2q2  +  p2q3  +  P3q2  +  p3q3 

The  non-negativity  of  the  entries  of  P,  and  the  marginality  conditions 
on  the  entries  of  P  coupled  with  P  being  PQD  yield  the  above  inequalities. 
These  inequalities  are  important  in  the  sense  that  if  there  are  four 
numbers  p22>  p23>  p32,  p33  satisfying  (1),  (2),  (3),  (4)  above,  then  we 
can  build  a  matrix  P  belonging  to  MpQp^pi’p2,p3’  ql’q2’q3^'  We  make 


precise  in  the  following  theorem. 


•4 
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Theorem  10  Let  p ^ »  p23»  P32*  p33  ^our  numbers  satisfying  the 
inequalities  (1),  (2),  (3),  (4)  above.  Define 

P21  =  P2  “  p22  “  p23 ’ 

P31  =  P3  "  P32  ”  P33* 

P^2  ~  ^2  —  P22  —  p32* 

P13  =  q3  “  P23  "  P33  * 

Then 

P1  "  P12  ~  p13  =  ql  P21  "  P31  =  pll»  say 
and  the  matrix  P  =  (p  )  e  »P2 »P3 5  ql »q2 ,q3)  * 

Proof.  The  proof  is  analogous  to  the  one  given  for  Theorem  7 . 
Determination  of  the  extreme  points  of  Mp  (p^^.p^;  qi,q2’q3^ 

As  explained  in  the  analysis  of  2  xn  tables,  inequality  (2)  gives 
rise  to  at  most  3  inequalities  equivalent  to  (1) .  Equivalent  inequalities 
are  given  below  for  each  of  the  cases. 


Case 

_A 

P3^ 

q3  i 

P2q3  +  p3q3 

and 

q3  ' 

P2  1  0 

(1) 

p3q3 

<  p33  <  p3a 

q3 

and 

(2) 

p2q3 

+  P3q3  1  P23 

+  ; 

p33  1 

q3 " 

Case 

_B 

P3  A 

q3  — 

p2q3  +  P3q3 

and 

0  <  c 

i3  -  p2  1  p3q3 

(1) 

P3q3 

-  p33  —  P3  ^ 

q3 

and 

(2) 

p2q3 

+  P3q3  i  p23 

+  ; 

P33  — 

q3 

Case 

_C 

p3  A 

q3  — 

P2q3  +  p3q3 

and 

P3q3 

<  q3  -  p2  <  p3  A  q3 

(l) 

P3q3 

-  P33  -  q3  " 

P2 

and 

(2) 

P2q3 

+  P3q3  1  P23 

+ 

P33  - 

P33  +  P2  > 

or 

(1) 

q3  ' 

P2  1  p33  -  P 

3  ^ 

q3  and 

(2)  p„q„  +  p,q-  <_  p,9  +  p,., 


<  q_. 
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Case  D 

P3  A 

q3  —  P2q3  +  P3q3  and  p3  A  ^3 

—  q3  ' 

(1) 

p3q3  -  P33  i  P3  A  q3  and 

(2) 

p2q3  +  p3q3  -  p23  +  P33  -  p2 

+  P33* 

Case  E 

P2q3 

+  p3q3  <  P3  /l  Q3  and  ^3  “  P2 

1  0 

(1) 

P3I3  i  P33  i  P2^3  +  p3q3  and 

(2) 

p2q3  +  p3q3  -  p23  +  p33  -  q3 

> 

(1)  p2q3  +  p3q3  <_  p33  £  p3  f\  q3  and 

(2)  p33  ±  p23  +  P33  ±  V 

Case  F  p2q3  +  P3^3  <  p3  A  ^3  anc^  ®  <  ^3  “  P2  ^3^3 
The  inequalities  here  are  the  same  as  those  of  Case  E. 


Case  G  p2q3  +  p3q3  <  P3  A  q3  an< *  p3^3  <  ^3  "  p2  <  p3  A  ^3* 
Assume,  without  loss  of  generality,  that  q3  -  p2  <  p2q3  +  p3q3' 


(1) 

P3q3 

1  p33  1  ^3  ~  P2  and 

(2) 

p2q3 

+  p3q3  1  p23  +  p33  <  p2  + 

or. 

(1) 

q3  ‘ 

p2  -  p33  1  p2^3  +  p3^3  and 

(2) 

p2q3 

+  p3q3  1  p23  +  p33  1  q3» 

or, 

(1) 

P2q3 

+  p3q3  i  P33  <_  P3  A  q3  and 

(2) 

P33  - 

-  P23  +  P33  -  q3 " 
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Case  H  p2q3  +  p3q3  <  P3  K  q3  and  P3  A  q3  <  q3  -  p2- 
(1)  p3q3  <_  p33  <_  p2q3  +  p3q3  and 

(^)  ^2^3  *^"  ^3^3  —  ?23  P33  —  ^2  ^33’ 

or, 

(1)  P2q3  +  P3q3  ±  P33  1  P3  A  q3  and 

(2)  P33  <_  p23  +  p33  1P2  +  P33- 

Under  each  of  the  Cases  A,  B  and  D,  the  inequality  (1)  virtually 
remains  the  same,  but  in  the  Cases  C,  E,  F  and  H,  the  inequality  (1) 
is  split  into  2  inequalities.  The  most  complicated  case  is  Case  G  in 
which  the  inequality  (3)  is  split  into  3  inequalities.  This  splitting 
of  (1)  influences  the  inequality  (2)  in  that  the  symbols  \/  and  f\ 
disappear  from  each  side  of  the  inequality  (2) . 

We  assume  that  we  have  Case  G  to  contend  with.  Note  that  we  have 

under  Case  G,  P393  <  ^3  ~  ^2  <  ^2q3  +  p3q3  <  p3  A  q3- 

Similarly,  the  inequality  (3)  gives  rise  to  at  most  3  inequalities 
equivalent  to  (1).  We  assume  that  here  also  we  have  the  most  complicated 
case  of  3  inequalities  equivalent  to  (1).  Assume,  without  loss  of 
generality,  that 

p3q3  <  P3q2  +  P3q3  <  P3  -  q2  <  P3  A  q3- 
Assume  also  that  P3q3  <  p2q2  +  p2q3  +  P3CI3  +  p3q3  <  p3  A  q3  ancl 

P3q3  <  q3  ~  P2  <  p3q2  +  p3q3  <  p3  "  q2  <  p2q3  +  P3q3 
<  p2q2  +  p2q3  +  p3q2  +  p3q3  <  p3  A  V 

If  some  other  order  prevails  among  these  seven  numbers,  we  can  analyze 
this  case  with  equal  facility.  We  are  dealing  with  the  maximum  number 
of  inequalities  that  can  arise  equivalent  to  (1)  induced  by  the 
inequalities  (2)  and  (3)  and  the  number  P2q2  +  ^2q3  +  p3q2  +  p3q3’ 


To  round  up  the  discussion  initiated  by  the  inequalities  (1),  (2) 
and  (3),  we  conclude  that  the  inequalities  (1),  (2)  and  (3)  are 
equivalent  to  the  following  sets  of  inequalities. 


[1] 

p3q3  i  p33  i  q3  -  P2» 

[ 1] (2) 

P2q3  +  P3q3  1  P23  +  P33  1  P2  +  P33.  and 

[1] (3) 

P3q2  +  P3q3  1  P32  +  P33  1  q2  +  P33> 

or. 

[2] 

q3  ‘  p2  —  p33  —  P3q2  +  p3q3  * 

[2] (2) 

p2q3  +  P3q3  i  P23  +  P33  1  q3»  and 

[2] (3) 

p3q2  +  p3q3  -  p32  +  p33  -  q2  +  p33‘ 

or. 

[3] 

p3q2  +  p3q3  1  p33  1  p3  '  q2‘ 

[3] (2) 

p2q3  p3q3  —  p23  p33  —  q3  ’ 

[3] (3) 

p33  -  p32  +  p33  -  q2  +  p33’ 

or. 

[4] 

p3  "  q2  -  p33  -  p2q3  +  p3q3  ’ 

[4] (2) 

p2q3  P3q3  —  P23  p33  —  q3  *  and 

[4] (3) 

p33  -  p32  +  p33  -  p3  * 

or. 

[5] 

p2q3  p3q3  —  p33  —  p2q2  p2q3  p3q2  p3q3* 

[5] (2) 

P33  ?22  p33  —  q3 »  and 

[5] (3) 

p33  -  p32  +  p33  -  p3  ’ 

or. 

[6] 

P2q2  +  p2q3  +  p3q2  +  p3q3  -  P33  ±  p3  A  q3  • 

[6](2) 

P33  1  P23  ■*"  P33  q3 »  an<i 

[6] (3) 

p33  -  p32  +  p33  -  p3  * 

We  observe  that  the  symbols  V  and  /\  have  disappeared  from  the 


inequalities  (2)  and  (3) . 
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Now,  we  have  to  fit  the  inequality  (4)  into  the  above  scheme  of 
inequalities  equivalent  to  (1) ,  (2)  and  (3) .  The  primary  goal  is  to 
get  rid  of  the  symbols  A  and  V  from  (4)  before  appending  it  to  each 
one  of  the  above  six  sets  of  inequalities.  Neutralization  of  the 
symbol  A  on  the  right  hand  side  of  the  inequality  (4)  can  be  effected 
in  the  same  way  as  its  done  in  the  inequalities  (2)  and  (3)  earlier. 
This  may  involve  further  splitting  up  the  inequalities  (2)  and  (3)  in 
each  and  every  set  of  inequalities  above.  Or,  A  can  be  retained  as  it 
is  in  (4) .  To  neutralize  the  symbol  ^  on  the  seft  hand  side  of  the 
inequality  (4) ,  we  introduce  a  new  inequality  intermediate  to  (2) ,  (3) 
on  one  side  and  (4)  on  the  other  side.  For  each  set  of  inequalities 
given  above,  one  can  determine  the  minimum  and  maximum  value  of 
p23  +  p32  +  p33*  Let  these  values  be  a£  and  a*  respectively.  There 
are  several  possibilities. 

Case  1  P2q2  +  p2q3  +  P3q2  +  p3q3  <  a*. 

In  this  case,  there  is  no  need  to  introduce  an  intermediate  inequality. 

Case  2  aj  <  p2q2  +  p2q3  +  p3^2  +  p3^3  <  a2 ' 

We  introduce  an  intermediate  inequality 

al  -  p23  +  p32  +  p33  -  p2q2  +  p2q3  +  p3q2  +  p3q3’ 
or, 

p2q2  +  p2q3  +  p3q2  +  p3q3  ±  p23  +  p32  +  P33  -  a2' 

Case  3  a£  <  p2q2  +  p2q3  +  p3q2  +  p3q3’ 

In  this  case,  there  is  no  need  to  introduce  an  intermediate  inequality. 
A  typical  set  of  inequalities  has  the  following  appearance. 

A  (1)  ax  <  p33  <  a2 

(2)  b1  <  p23  +  p33  <  b2 

(3)  c1  <_  p32  +  p33  £  c2 

(4)  dx  <  p22  +  p32  +  p23  +  p33  1  d2; 


e1  <_  P22  +  p23  +  P32  +  p33  £  e2. 

In  the  first  lot  above,  b^,  b2,  c^  and  c2  could  be  either  constants  or 
linear  functions  of  p33»  The  same  remarks  apply  to  b^ ,  b2>  c^  and  c2 

in  the  second  lot.  In  the  first  lot,  d^  =  p23  +  p32  +  p33  or 

P2q2  +  V2*2  +  P3q2  +  p3q3  and  d2  =  P2  +  p32  +  P33  or  q2  +  p23  +  p33. 

In  the  second  lot,  d^  =  minimum  value  of  p23  +  p32  +  p33  subject  to 

(1),  (2)  and  (3)  or  p2q2  +  p2q3  +  p3q2  +  p^  and  d2  =  p2q2  +  p2q3  +  p3q2  +  p3q 

or  the  maximum  value  of  p23  +  p32  +  p33  subject  to  (1),  (2)  and  (3). 

Further,  e;l  =  p23  +  p32  +  p33  or  p2q2  +  p2q3  +  p3q2  +  p3q3  and 

e2  =  P2  +  P32  +  P33  °r  q2  +  p23  +  p33  * 

In  conclusion,  we  have  some  sets  of  inequalities  where  each  set  is 

one  of  the  two  types  A  or  B  described  above.  These  sets  of  inequalities 

are  jointly  equivalent  to  (1),  (2),  (3)  and  (4). 
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In  order  to  find  the  extreme  points  of  HpQD(p1  ,p2>p3;  q-^q^q-j), 

we  look  at  the  following  problem.  Let  A  be  the  collection  of  all 
vectors  (p22,  p23,  p32,  p33)  satisfying 
(1)  ax  <_  p33  <_  a2, 

(2>  bi  i  P23  +  p33  -  b2* 

(3)  cx  p32  +  p33  _<  c2, 

(4)  dx  £  p23  +  p32  +  p33  <_  d2  and 

(5)  p22  +  p23  P32  P33  —  e2  * 

where  a^  and  a2  are  given  fixed  numbers;  b^,  b2,  c^,  c2  are  functions 

of  p33;  d^,  d2  are  fixed  numbers;  and  e^,  e2  are  functions  of  P23 >  P32 

and  P^3*  The  set  A  is,  obviously  a  compact  convex  set  with  a  finite 

number  of  extreme  points.  We  now  describe  a  method  of  generating 

extreme  points  of  the  set  A.  Let  F»  »  .  .  .  denote  the  collection 

61*62*°3,64»°5 

of  all  vectors  (p22»  P23»  P22>  P33)  *n  ^  for  which 

equality  holds  in  (i)  on  the  left  of  (i)  if  <5^  =  1, 

equality  holds  in  (i)  on  the  right  of  (i)  if  6  *>  -1, 

strictly  inequality  holds  in  (i)  if  6^  =  0, 

for  i  =  1,2, 3, 4, 5.  Each  6^  takes  one  of  the  values  1,  -1,  0.  There 

3“*  =  243  pairwise  disjoint  sets  F.  *  *  «  »  's  whose  union  is  A. 

1  2  3  4  5 

Some  of  these  sets  could  be  vacuous.  Assume  that  for  each  6^,62,<$3,<5^,S3 
there  exists'  <$,  such  that  Ff  .  .  r  .  ^  0.  This  assumption  is 


=  +1  or  -1 


satisfied  in  the  original  problem  under  investigation.  We  define  the 
dimension  of  each  of  these  sets  as  follows.  We  associate  the  following 
vectors  with  each  of  the  five  inequalities  (1),  (2),  (3),  (4)  and  (5). 


i|  iji  my  n  ijij  ij  r-  r. ' 
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rj  ■  y  i  j  *  : »  i  " 


Inequality 

(1) 

(2) 

(3) 

(4) 

(5) 


Vector 
(0,0, 0,1) 
(0,1, 0,1) 
(0,0, 1,1) 
(0,1, 1,1) 
(1,1, 1,D 


Let  6- ,<5_ ,6_ ,6.  and  5-  be  fixed.  The  dimension  of  F.  ....  is 
1  2*  3*  4  5  fi1,<S2,63,<S4,<S5 

defined  to  be  the  number  of  4  -(the  number  of  linearly  independent 

vectors  among  (i) 's  above  for  which  6^  ■  +1  or  -1).  For  example,  the 

set  F.  .  .  i  has  dimension  0.  The  F1  -  Tin  ^as  dimension  1. 


Each  of  the  sets  F 


61’62’63*64’65 


is  compact  and  convex.  These 


sets  have  the  following  property. 

Theorem  11  Every  member  of  any  set  of  dimension  i  is  a  convex  combina¬ 
tion  of  members  of  some  sets  of  lower  dimension  for  i  «  1,2, 3, 4. 

Proof  The  proof  in  all  its  details  is  very  lengthy.  We  merely  sketch 
the  details.  Take  any  set  of  dimension  1,  Fn  .  .  .  n,  say.  Look  at 

the  sets  F..  .  1  .  1  and  F.  .  .  ,  . .  These  two  sets  are  singleton 
XfX,— X,— X,— X  XjXj^X,— x,x 

sets.  If  (p22’P23,p32’p33)  6  Fl,  1,-1, -1,0.’  (<*22’q23’q32’q33)  e 

Fl, 1,-1, -1,-1  and  (s22’s23’s32’s33)  E  Fl, 1,-1, -1,1’  then  p33 =  q33 =  s3 

P23  m  *^23  m  s23’  p32  *  q32  *  832  ^**d  q22  ^  p22  ^  822*  can  write 
(p22,p23’p32,p33^  as  a  convex  combination  of  (q22 ,q23 ,q32 ,q33^  and 
^S22*823’S32’S33^ *  Similarly,  every  member  of  a  set  of  dimension  2 
can  be  written  as  a  convex  combination  of  some  two  members  of  a  well- 
chosen  set  of  dimension  1.  The  same  argument  can  be  carried  out  for 
sets  of  higher  dimension. 

As  a  consequence  of  the  above  result,  we  can  identify  the  extreme 


points  of  A. 


Corollary  12  The  set  of  extreme  points  of  A  is  contained  in  the  union 


of  all  sets  of  dimension  0. 

Theorem  11  and  Corollary  12  are  useful  in  enumerating  extreme 
points  of  MpQD(Pi»P2*P3»qi»q2»q3^ *  For  each  set  of  inequalities,  using 
Theorems  11  and  12,  enumerate  the  extreme  points  of  the  corresponding 
set  A  defined  as  above.  The  totality  of  all  these  extreme  points  pooled 
together  from  each  set  of  inequalities  constitutes  the  set  of  all 
extreme  points  of  MpQD(p1,p2>p3;q1,q2,q3) . 

We  illustrate  the  above  procedure  with  some  examples. 

Example  13 

Let  p^  =  p2  =  p3  *  q^  =  q2  =  q3  =  1/3 
Determining  inequalities  are: 


(1) 

p3q3 

—  p33 

1  p3  A  q3» 

(2) 

p33 

V (P293 

+  p3q3)  1  p23 

+  P33  i  <*3 

A 

(p2 

+  P33)  , 

(3) 

P33 

V  (p3q2 

+  p3q3)  <  p32 

+  P33  i  P3 

r\ 

(q2 

+  P33), 

(4) 

(p23 

+  p32 

+  p33)  v  (p2q2 

+  P2q3  +  P; 

3q2 

+ 

P3q3) 

1  P22  +  p23  +  1 

P32  +  P33 

i  (p2  +  p32  +  P33)  ^  Cq2  +  P23  +  ' 


After  substitution,  these  inequalities  become 

...  1  1 

^  9  -  p33  —  3’ 

(2)  p33  v/  (-|)  <_  p23  +  P33  1  J  A  (j  +  P33)  » 

(3)  P33  \/  (-9)  1  P32  +  P33  -  3  ^  ^3  +  p33^  ’ 

(4)  (p23  +  p32  +  p33)  V (f)  1  p22  +  p23  +  p32  +  p33 

-  %  +  p32  +  p33^  ^  ^ 3  +  p23  +  P33 

The  inequalities  (2)  and  (3)  induce  a  partition  of  the  inequality  (1) 


). 


The  inequalities  (1),  (2)  and  (3)  are  equivalent  to 

U1U>  |iP33if 

[i]  (2)  jiPjj  +  P33  i  3 

[1]  (3)  f  i  P32  +  P33  1 
and 

(1)  y  1  P33  1  y 

[2]  (2)  p33  <_  p23  +  P33  i  y 

[2]  (3)  p33  <  p32  +  P33  1  -J 

The  maximum  value  of  p23  +  p32  +  p33  when  p23>  p32 ,  p33  satisfy 
C 1 3 (1 ) ,  [ 1 ] (2 )  and  [1](3)  is  y.  The  maximum  value  of  p23  +  p32  +  p33 
when  p23<  p32,  p33  satisfy  [2 ] (1) ,  [2 ] (2)  and  [2](3)  is  4.  Inequality  (4) 
can  be  appended  to  each  one  of  the  sets  above  giving  rise  to  3  sets  of 
inequalities  equivalent  to  (1) ,  (2) ,  (3)  end  (4) . 

[1]  (1)  P33  if 

[!]  (2)  y  <_  p23  +  p33  1  y 

[1)0)  |  1  P23  +  P33  1  y 

[11  (4)  y  £  P23  ’’’  P32  P33  —  "9 

[1] (5)  g  £  P22  P23  ^  P32  P33  —  ^  P32  P33)  ^  (3  P23  P33)  » 

or, 

I21(1)  i  i  P33  -  I 

[2]  (2)  y  1  P23  +  P33  I  2 

12H3)  |1P32  +  P33<I 

[2!f4!  ~9  -  P 21  +  P 32  +  P33  -I 

12  1,1  Pj  (  +  P  i 2  +  P33  1  P22  +  P2 3  +  P32  +  P33 

1  ^  ~T  +  Pot  P  o  o )  ^  ^  t  Poo  ^  P  o  t  ^  » 
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I3](l)  f  1  p33  ly 

[3]  (2)  p33  <_  p23  +  p33  <_  y 


[3] (3) 
f 3] (4) 


P33  -  p32  +  p33  -  i 

■9  -  p22  +  p23  +  p32  +  p33  -  %  +  p23  +  p33^  ^  ( T  +  P32  + 


p33 


Now,  in  each  of  the  sets  [1]  and  [2]  and  in  each  of  the  inequalities 
(1) ,  (2),  (3)  and  (5),  we  set  the  central  expression  either  equal  to 
the  quantity  on  the  right  or  on  the  left  and  then  solve  the  resultant 
set  of  equations  in  the  unknowns  P22»  P23*  p32  anc*  p33  ma^-^n8  sure 
that  the  inequality  (4)  is  satisfied.  Using  these  numbers,  we  build 
a  matrix  belonging  to  Mpqj^yj  y»  y’>  y»  y»  y)  •  Avoiding  repititions,  we 
obtain  the  following  matrices. 

1 
9 

1 
9 

1 
9 


<-H|cr>  o  cn|(Ti  O  i— i  |  cr\  cmI'J'  O  i-h|<3'.  csijCT' 


rz 


9 


In  the  third  lot  of  inequalities,  set  the  central  expression  in 
each  inequality  either  equal  to  the  quantity  on  the  right  or  on  the 
left  and  then  solve  the  resultant  equations  in  $22’  P23»  P32  and  P33* 
We  obtain  the  following  matrices  after  omitting  repititions. 


1  2  1 

r  3  1 

9  9  0 

18. 

0 

0 

t|  O' 

2  1. 

„  3 

9  9  0 

0  y  0 

3 

3 

0  0  4 

O 

O 

9 

9  j 

Note  that  the  matrix  14  is  a  convex  combination  of  the  matrices 

16  and  18.  All  these  seventeen  matrices  with  the  omission  of  matrix  14 

..,..,111111, 

are  the  extreme  points  of  M^Cy,  y,  y  y,  y,  y) • 

It  is  interesting  to  compare  the  extreme  points  of  the  set 

MpQI)(y,  f’  V  T»  y*  wlth  the  set  M(y>  y>  J»  y»  y«  y)  of  all  matrices 

of  order  3x3  with  nonnegative  entries,  row  sums  y,  y,  y  and  column 

sums  y,  y,  y.  The  set  M(y,  y,  y;  y,  y,  y)  has  only  6  extreme  points 

given  by 


Bi 

■ 

■  I 

I 

■ 

H 

.  8 

1 

■9 

H  ■■  1 

:  88 

HI 

mV 

■  ' 

H^CB 

■  ■ 

■ 

|  | 

^  X 

9  C";i 

■  ■ 

B 

■ 

II 

:  _  8 

SB 

1 

■  I 

n 

1 

■  ■ 

8 

■  8 

8 

■Mm 

B 

91 

■ 

1 

1 

8 

B 

B 

1 

ml 

1 

|B»j 

■■ 

Mm \ 

Note  that  the  set  y.  "3!  "J»  y»  y)  is  a  subset  of 

M("3 »  y,  3-;  -J.  -3,  j)  and  only  one  extreme  point  of  M(-j,  j,  j;  j,  j,  -) 
Is  an  extreme  point  of  J,  y  J»  3>  • 

Example  14 

.  .  _  1  _  1  _  1  _  1  _  1  ,  3 

Let  2  9  g>  ^2  8  ^2 

The  determining  inequalities  are: 

(1)  p3q3  <  p33  <  p3  A  13 

(2)  P33  ^  ^2*^3  p3q3^  —  p23  P33  —  q3  A  ^2  p33^ 

(3)  P33  V  ^ P 3*^-2  p3q3^  —  P32  ^33  —  p3  A  ^2  p33^ 

(4)  (p23  +  p32  +  p33)  V  (p2q2  +  p2q3  +  p3q2  +  p3q3) 

-  p22  +  p23  +  P32  +  P33  —  (p2  +  p32  +  p33^  A  (q2  + 
After  substitution,  these  inequalities  become 

(1)  re  1  p33  i  7 

(2)  p33  V  (  16  ^  -  P23  +  p33  -  T  ^  +  p33^ 

(3)  p33  \/  (  32  )  i  P32  +  P33  1  "4  A  (g-  +  P33) 

2 1 

(4)  (p23  +  p32  +  p33)  V  (-^r)  £  p22  +  p23  +  p32  +  p33 

-  (' 2  +  p32  +  p33^  ^  %  +  p23  +  P33^  * 

The  inequalities  (1),  (2)  and  (3)  are  equivalent  to  the  following  two 


sets  of  inequalities. 


t2]  O)  P33  i  P32  +  P33  i  4 


The  maximum  and  minimum  values  of  p23  +  P32  +  P33  are  tabulated 


below  for  each  of  the  above  two  sets. 

Sets  of  inequalities  Expression 


[1] 

P23  +  p32  +  P33 

[2] 

p23  +  p32  +  p33 

Max imum  Minimum 

24  18 


We  append  the  inequality  (4)  to  each  one  of  the  above  sets  after 

18  21 

splitting  the  range  of  p23  +  p32  +  p33  into  two  parts  and 

21  24 

[~32*  ^2”] .  We  obtain  four  sets  of  inequalities  equivalent  to 
(1),  (2),  (3)  and  (4). 

UKl)  32  -  P33  <32 

[i] (2)  ^2  -  p23  +  P33  -  32  +  p33 

[i]  (3)  P32  +  P33  <.32 

mz/N  18  .  .  21 

[  ](  )  32  —  p23  +  P32  +  P33  -  32 

UK5)  32  —  p22  +  p23  +  p32  +  P33  -  (l2  +  P32  +  p33^  ^  ^32  +  p 


[2]  (1)  1  P33  1  32 

[2]  (2)  -J2  1  P23  +  p33  -  ^2  +  p33 

[2]  (3)  ^  <  P32  +  P33  1  -32 

[2] (4)  32  i  P23  +  P32  +  p33  -  32 

[2K5)  P23  +  P32  +  P33  1  ?22  +  p23  +  p32  +  p33 


-  (32  +  p32  +  P33^  A  (32  +  p23  +  P33^ ’ 


[3] (5) 


32  -  P22  +  p23  +  P32  +  P33  -  (32  +  p32  +  p33^  ^  (32  +  p23  +  p335  ’ 


[4](1) 

*. 

[4] (2) 
[4] (3) 
[4] (4) 
[4] (5) 


32  -  P33  -  32 

18  ^  16 

32  -  P23  +  P32  -  32  P33 

g 

p33  -  P32  +  p33  -  32 

21  24 

32  —  p23  +  pi?  +  pv*  —  n 


32 


33  -  32 


23 


+  P32  +  Pri  1  P??  +  P?i  +  P-i?  +  P 


33 


22 


23 


32 


J33 


-  (32  +  P32  +  P33^  ^  (32  +  p23  +  P33) * 


In  each  set,  for  each  of  the  inequalities  (1),  (2),  (3)  and  (5), 
we  set  the  central  expression  equal  to  either  the  quantity  on  the  left 
or  on  the  right  and  then  solve  the  resultant  equations  in  P22»  P23>  P32 

and  after  making  sure  that  the  inequality  (4)  is  satisfied.  Using 
p22 >  P23»  P32»P33>  one  can  build  p  in  ql’q2,q3)  in  a 

natural  way.  These  matrices  are  given  below. 


5 


1 

1 

6 

32 

32 

32 

3 

3 

10 

32 

32 

32 

0 

0 

8 

32 

2 

0 

6 

32 

32 

2 

4 

10 

32 

32 

32 

0 

0 

8 

32 

These  are  the  extreme  points  o 


V 


4.  On  the  power  of  tests  for  testing  independence  against  positive 


quadrant  dependence 

In  this  section,  we  apply  the  results  of  the  previous  sections  to 
study  power  of  tests  in  the  context  of  contingency  tables.  The 
problem,  basically,  is  the  following.  Let  X  and  Y  be  two  random  • 
variables  and  each  takes  finitely  many  numerical  values.  The  distribu¬ 
tions  of  each  of  X  and  Y  are  known  but  their  joint  distribution  is 
unknown.  For  simplicity,  assume  that  X  takes  values  l,2,...,m  with 
probabilities  Pi>p2»****Pm  respectively  and  Y  takes  values  l,2,...n 
with  probabilities  respectively.  Let  p^  =  P(x=i,  y  =  j}, 

i  =  1,2, ... ,m  and  j  =  1,2, ... ,n.  We  say  that  X  and  Y  are  strictly 
positive  quadrant  dependent  if  X  and  Y  are  positively  quadrant  dependent 
but  X  and  Y  are  not  independent .  We  want  to  test  the  validity  of  the 
null  hypothesis 


against 


Hq:  X  and  Y  are  independent 


H^:  X  and  Y  are  strictly  positive  quadrant 


dependent 


based  on  a  random  sample  of  size  N  on  (X,Y) .  Let  n^  =  number  in  the 
sample  with  x  =  i  and  y  =  j,  i  =  1,2, ...,m  and  j  =  1 , 2 , . . .  ,n .  The 
data  can  be  tabulated vin  the  form  of  a  contingency  table. 


\  y 

x  \ 

1 

2 

•  «  • 

n 

1 

nll 

n12 

•  •  • 

ni 

In 

2 

n21 

n22 

•  •  • 

c 

fsl 

c 

• 

• 

• 

; 

•  •  • 

• 

m 

• 

nml 

nm2 

•  *  • 

• 

n 

mn 

m  n 

Here,  J  I  n  =  N. 
i=l  j-1  J 


Let  T  be  any  test  proposed  for  testing  against  for  a  given 
level  of  significance,  a.  T  is  a  function  of  the  random  variables 
n^  ,  i  =  l,2,...,m;  j  =  l,2,...n.  The  performance  of  the  test  T  can  be 
judged  based  on  its  power  function.  More  precisely,  let 
MpQD(Pi»P2»  •  •  •  »Pm5  '  *  *  ,c*n^  t*ie  collection  of  all  matrices  P  of 

order  m  *  n  with  non-negative  entries,  row  marginal  sums  Pi>P2 » • • • »Pm* 
column  marginal  sums  q^,q^,...,qn  and  positive  quadrant  dependent. 

Let  Pq  be  the  matrix  of  order  m  x  n  whose  ( i , j ) t^1  element  is  given  by 
p^q^  ,  i  =  1,2,..., m  and  j  =  l,2,...,n.  The  parameter  space  of  the 
above  problem  is  MpQD(p1  ,p2 , . . .  .p^;  q^^  ,q2  , . . .  ,q^)  .  The  hypotheses 
Hq  and  translate  into  as  follows.  Let  P  stand  for  the  generic 

symbol  in  mPqD(p1»P2 . Pm5  *  ••  •  > cln)  deemed  as  the  joint 

distribution  of  X  and  Y.  The  relevant  hypotheses  are: 

Ho:  p  =  po 
H1:  P  *  p0 

The  power  function  $(.)  of  T  is  defined  by 

Bt(P)  =  Pr  {T  rejects  Hq  /  the  joint  distribution  of 

X  and  Y  is  P} , 

P  e  MPQD(pl,P2’"  ,,Pm;  ql*q2* ’ *  *  ,qn^  ’ 

Obviously,  &t(Pq)  =  a.  The  following  results  advocate  the  necessity 
of  enumerating  the  extreme  points  of  >P2 * • • • »Pm»  qi ,q2 * ’  '  * ,qn^ * 

Theorem  4.1  Let  P  e  MpqD^Pi»P2» • • • >Pm»  ql ,q2 * ' ’ ,qn^  be  3  convex 
combination  of  some  matrices  P^  ,P2 ,  . . .  ,P^  in  ^pq^P  p  *P2  »  •  •  •  >Pra>  » CI2  »  • 


p  =  l  «i  Pi 
i=l 

k 

for  some  a.,a„,...,a,  0  with  £  a  =  1.  Then 

l  z  x  i=l  1 


8t(P)  =  l  <*<  B  (Pi) 
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Proof.  Obvious. 

The  above  theorem  establishes  the  fact  that  £(.)  is  an  af f ine 

function  on  the  compact  convex  set  MpQD(P1»P2 . Pm5  qx»q2 . • 

Corollary  4.2  Let  E  (p1,P2> • . . ,Pm;  q1,q2,...,q  )  =  E,  say,  be  the  set  of 
all  extreme  points  of  mPqD<P1»P2> *  * • »Pm;  q1,q2> • • • ,qn) .  Let  E 
consist  of  some  k  matrices  Pi»p2 ’ * ’ * ,Pk*  Then  for  any  P  in 
MpQD(Pi»P2» • • • »Pm5  qi’^2’ *  * ' ,qn^ »  is  a  convex  combination  of 

eT(p1),  eT(p2),...,  eT(pk). 


The  above  corollary  provides  us  a  criterion  by  which  we  can 
compare  the  performance  of  any  two  given  tests  T^  and  T2  proposed  to 
discriminate  between  Hq  and  H^.  If  P^,P2>...,P^  are  the  extreme 
points  of  MpqD(Pi>P2»  •  •  •  >Pm»  ^l,^2»  *  *  •(ln)  »  then  we  compare  the  two 
sets  of  numbers 


{  0T1(P1)  ’  BT1(P2)  »*••*  6Ti(Pk)  1 
and 


*  6T2^P1^  *  6T2^P2^  **••»  ST2^Pk^  * 

for  a  showdown  between  T^  and  T2  as  to  which  one  is  more  effective  in 
discriminating  between  Hq  and  H^. 

In  a  paper  by  Bhaskara  Rao,  Krishnaiah  and  Subramanyam  [6],  a  detailed 
comparison  of  performance  of  several  tests  is  carried  out  for  some  cases 
of  contingency  tables. 


.\ 


5.  Comparison  of  tests 

In  this  section,  we  actually  embark  on  comparing  the  performance 
of  several  tests  for  testing  the  null  hypothesis  Hq  against  the  alternative 
specified  in  the  previous  section  in  the  case  of  2x2  and  2x3 
contingency  tables.  First,  we  deal  with  the  2x2  case. 


We  now  discuss  the  power  function  of  the  test  based  on  Goodman- 
Kruskal ’ s  Gamma  Ratio.  Let 


be  the  joint  distribution  of  X  and  Y  with  specified  marginal  distributions 
P1'P2  and  qi'q2  resPectiv®ly-  The  Goodman-Kruskal’ s  Gamma  Ratio,  y, 
is  defined  by  (see  Agresti  [l ] ) 


P11P22 


P12P21 


P11P22  +  P12P21 


The  following  properties  of  -y  are  easy  to  check. 

(i)  --1  <  y  <  1. 

(ii)  X  and  Y  are  independent  if  and  only  if  y  =  0. 

(iii)  X  and  Y  are  positive  quadrant  dependent  if  a'nd  only  if  y  >_  0. 

(iv)  X  and  Y  are  strictly  positive  quadrant  dependent  if  and  only 

if  y  >  0. 

From  the  above  properties,  we  can  assert  that  y  is  a  measure  of  positive 
quadrant  dependence  between  X  and  Y.  (For  the  validity  of  (iii)  ,  we  have 
used  the  fact  that  the  marginal  distributions  are  specified. )  We  also 


observe  the  following.  The  extreme  points  of  the  convex  set  Mpqd^pi ,p2;C5l ,q2^ 


are  given  by 


Plqi  Piq2 


P2qi  P2q2 


assuming  that  P2  ^  q2  =  P2* 
of  X  and  Y  is  P  ,  y  =0 


and  P2  = 

Under  P  ,  i.e.,  the  joint  distribution 
and  under  p2 '  ^  ~  1  - 


qi  q2~P2 


(5.1) 


A  natural  estimator  of  y  is  given  by 


Y  = 


nlln22  ~  nl2n21 
nlin22  +  ni2n21 


(5.2) 


Test  of  H  against  H,  based  on  y  :  T. 

- - o  — -  - 1_ - — 

Reject  in  favour  of  if  and  only  if  y  >  a, 

f 

where  a  is  chosen  satisfying  Pr  (y  >_  a  /  H^)  =  a»  the  given  level 

of  significance.  One  can  work  out  the  power  function  of  this  test  explicitly 
Let 

B_(P)  *  Pr(the  test  T  rejects  H  /  P) 

T  o 

*  Pr  (y  a  /  P) 

for  P  in  Mpqd^P^  'P2'ql'q2^  '  w^ere  the  above  probability  is  computed 
when  X  and  Y  have  the  joint  distribution  P.  As  has  been  pointed  out  in 
Section  4,  6,p(P)  is  connected  to  gT(P^)  and  6t^p2^  as  follows.  Write 

P  =  ip1  +  (1-1)P2 

for  some  O  <_  \  <_  1.  Then 

eT(P)  =  ABtTP1)  "+  I1=i)6T(P2)' 

Obviously,  $  (P  )  =  o.  Under  P  ,  n  n  =  0  almost  surely  and 
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consequently  Y  -  1  almost  surely.  This  implies  that  B^(p^)  *  Thus 
we  have  proved  the  following  theorem. 

Theorem  5.1  The  power  function  S  (  )  of  the  test  T  of  size  a  is 
given  by 

6  (P)  *  Xa  +  (i-X) 

T 

for  P  in  MpgD(P1'P2;qi'q2)  '  where  p  =  X'pl  +  U--X)P2- 

The  above  result  yields  the  following  result. 

Theorem  5 . 2  Let  T*  be  any  test  proposed  for  testing  Hq  against 
specified  in  the  previous  section.  Then 

BT<P>  i  BT.<P> 

for  every  P  in  MPQD (PX -P2 ;qi 'q2} ' 

Thus  we  see  that  in  the  case  of  2x2  contingency  tables,  the 
test  based  on  Goodman-Kruskal ’ s  Gamma  Ratio  is  superior  to  any  test  one 
can  propose  for  testing  the  hypothesis  of  independence  against  the  hypothesis 
of  ;strict  positive  quadrant  dependence. 

Let  us  examine  how  the  test  based  on  Goodman-Kruskal 1 s  Gamma 
Ratio  fares  against  another  test  proposed  in  the  case  of  2x3  contingency 
tables. 

Let  the  joint  distribution  of  X  and  Y  be  given  by 


*1! 

P12 

P13 

P1 

P21 

P22 

P23 

P2 

qi 

q2 

q3 

1 

v* 


where  p^  ,p^  and  ajre  t^le  niarginal  distributions  of  X  and  Y 

respectively.  The  Goodman-Kruskal* s  Gamma  Ratio,  y  ,  is  defined  by 


where  II 


C  *  P11(P22  +  P23^  +  P12P  23  “d 


and  n 


P13(P21  +  P22}  +  P12P21 ‘ 


See  Agresti  £l,  p.l6o].  We  propose  another  measure  of  association  between 
X  and  Y  as  follows. 


(P11P22  '  P12P215 


(P11P23  '  P13P21) 


(p12P23  P13P225 


PlP2qiq2 


P:P2qiq3 


PlP2q2q3 


of 


This  measure  of  association  is  motivated  by  the  definition^Spearman1 s  rho 


for  2x2  contingency  tables.  The  following  properties  of  y  and  < 
are  easy  to  verify. 


1. 

-1  <_  y  <_  1. 

2. 

If  X  and  Y 

are  independent,  then 

Y  =  0. 

3. 

If  X  and  Y 

are  positive  quadrant 

dependent,  then 

4. 

0  <,  K  <_  1. 

5. 

k  =  0  •  if  and 

only  if  X  and  Y  are 

independent. 

One  can  build  tests  based  on  the  Gamma  Ratio  and  the  kappa 


criterion. 


Test  based  on  the  Gamma  Ratio  :  T 


Let  C  *  the  total  number  of  concordant  pairs  in  the  data 
*  nll(n22  +  n23)  *  n12n23 

and  D  =  the  total  number  of  discordant  pairs  in  the  data 
=  nl3(n2i  +  n22)  +  n12nn> 

An  estimate  of  y  is  given  by 

*  _  C  -  D 

Y  v~  C  +  D 

Test  based  on  y  is  given  by  : 

Reject  the  null  hypothesis  in  favour  of  if  y  >  a, 

where  a  is  given  by  Pr  (y  >_  a  /  Hq)  =  o»  the  given  level  of  significance 
Let  6T(  )  be  the  power  function  of  the  test  T. 


Test  based  on  the  Kappa  value  :  T^_ 


An  estimate  of  <  based  on  the  given  data  is  given  by 


(nlin22  ~  n!2n21)  +  (nHn23  ~  n!2n21} 


(n!2n23  '  ni3n22 


N  Pi*2qiq2 


n  p1p2q1q3 


N  P,P,q9q, 


Test  based  on  <  is  given  by  : 

A 

Reject  in  favour  of  if  and  only  if  *c  >_  a, 

where  a  is  given  by  Pr(<  >_  a  /  H  )  =  a>  the  given  level  of  significance 


Let  (  )  be  the  power  function  of  the  test  T*. 


